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ABSTRACT 

In attempts to explain dark energy, a number of models have been proposed in which 
the formation of large-scale structure depends on the local environment. These models 
are highly non-linear and difficult to analyse analytically. A^-body simulations have 
therefore been used to study their non-linear evolution. Here we extend excursion 
set theory to incorporate environmental effects on structure formation. We apply the 
method to a chameleon model and calculate observables such as the non-linear mass 
function at various redshifts. The method can be generalized to study other obervables 
and other models of environmentally dependent interactions. The analytic methods 
described here should prove useful in delineating which models deserve more detailed 
study with JV-body simulations. 
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1 ITNTRODUCTION 

One of the most challenging questions in contemporary 
physics is the nature of the dark energy, which is believed 
to be driving the accelerating expansion of the Universe 
(Riess et al. 1998; Perlmutter et al. 1999). Copeland et al. 
(2006) present a comprehensive review of theoretical models 
to explain the apparent acceleration of the Universe. How- 
ever, at present there is no compelling evidence for any new 
physics other than the addition of a cosmological constant 
to the Einstein field equations. 

Models of dark energy can be broadly placed into two 
categories. In the first, the dark energy affects the expan- 
sion rate of the Universe but does not interact directly 
with the dark matter. Examples of this type of model in- 
clude the standard ACDM paradigm and quintessence mod- 
els (Wang et al. 2000). In the second category, the dark en- 
ergy and matter (both dark and baryonic) interact with each 
other with an interaction strength which may depend on the 
local environment. Examples include the chameleon coupled 
scalar field model (Khoury & Weltman 2004; Mota & Shaw 
2007), f(R) gravity (Carroll et al. 2005), the environmen- 
tally dependent dilaton model (Brax et al. 2010) and also 
the symmetron model (Hinterbichler & Khoury 2010). 

The possibility of environmentally dependent interac- 
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tions needs to be considered when relating laboratory mea- 
surements to cosmological scales. Consider, for example, a 
scalar field coupled to matter. The scalar field could mediate 
a 'fifth force' between matter particles. Current laboratory 
experiments and Solar System tests have shown that such 
a fifth force must be either extremely weak or short-range 
(less than about a millimetre) (Will 2006). However, it is 
possible that the strength and range of the fifth force de- 
pend on the environment so that locally, where the matter 
density is high, it is strongly suppressed, and it is restored 
in empty environments. In this situation, laboratory exper- 
iments cannot constrain a fifth force that may have obser- 
vational consequences on cosmological scales. 

Analytical models of structure formation on galaxy and 
cluster scales are notoriously difficult even in the case of 
standard Newtonian gravity. The evolution of structure in 
models with environmentally dependent interactions is even 
more complicated because the fifth force itself is highly 
non-linear. Consequently, studies so far have relied on full 
A-body simulations (Li & Zhao 2009, 2010; Li & Barrow 
2011; Li 2011; Li, Mota & Barrow 2011; Li et al. 2011; 
Zhao, Li & Koyama 2011a; Brax et al. 2011; Davis et al. 
2011; Oyaizu 2008; Oyaizu et al. 2008; Schmidt et al. 2009). 

However, large A-body simulations require supercom- 
puting resources and are time consuming. They can be jus- 
tified for testing physically well motivated models such as 
the ACDM model, which contains few parameters, many 
of which are now well constrained experimentally, see e.g., 
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Komatsu et al. (2011). Models with a fifth force, on the 
other hand, span a wide parameter space reflecting our lack 
of knowledge of the underlying physics. It is difficult to sam- 
ple a large parameter space using full iV-body simulations, 
hence the need for an analytic description of structure for- 
mation that can, at least, isolate regions of parameter space 
that merit further investigation using simulations. 

Semi-analytical models, such as excursion set theory 
(see Zentner (2007) for a recent review), have been devel- 
oped as alternatives to full iV-body simulations and shown 
to agree with the latter well. The excursion set approach has 
been generalised to some non-standard structure-formation 
scenarios (Martino et al. 2009; Parfrey et al. 2011). How- 
ever, these studies do not consider the case of environmen- 
tally dependent interactions. 

The aim of this paper is to generalise the excursion 
set approach to take account of environmentally dependent 
interactions explicitly. As we will see, non-linear collapse 
of structures could be very different in different environ- 
ments, and indeed the environments themselves evolve in 
time as well. We will first specify the environments using 
what we call the fixed-scale approximation, then use a sim- 
plified model to study spherical collapse within these envi- 
ronments. We then calculate observable properties by aver- 
aging over the distribution of environments. In this paper 
we have chosen the chameleon model as a working example, 
but the methods developed are more general and with suit- 
able changes can be applied to other models with environ- 
mentally dependent interactions. The theoretical framework 
developed here can therefore be used to quickly estimate the 
parameter ranges of any specific theory that may have inter- 
esting (and potentially testable) consequences on structure 
formation. 

The layout of this paper is as follows. We introduce the 
basic formulae for a chameleon-like coupled scalar field (our 
working example) in § 2 and summarise the spherically sym- 
metric solutions which will be used later to study the spher- 
ical collapse of overdensities. § 3 presents the main results of 
this paper. We introduce the traditional excursion set theory 
in § 3.1 and in § 3.2 we show how the environmental depen- 
dence in the chameleon model can be approximated using 
only two variables. § 3.4 describes a generalised spherical 
collapse model in which an overdensity collapses inside an 
evolving environment. Finally in § 4 we make an application 
of the generalized excursion method to a range of chameleon 
models. Our conclusions are summarized in S 5. 



are presented here for completeness and to establish the no- 
tation used in later sections. 

We start from a Lagrangian density 



C 
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in which R is the Ricci scalar, k = SirG with G being the 
gravitational constant, £dm and Cs are respectively the 
Lagrangian densities for dark matter and standard model 
fields, ip is the scalar field and V(ip) its potential; the cou- 
pling function G'(ip) characterises the coupling between ip 
and dark matter. Given the functional forms for V(ip) and 
C(ip) a coupled scalar field model is then fully specified. 

Varying the total action with respect to the metric g a b, 
we obtain the following expression for the total energy mo- 
mentum tensor in this model: 

iv c v^-v%>) 



T ab = VaV3V(,</3 - g ab 



+ C{<p)T° b M + Ta b , 



where T^ M and T^ b are the energy momentum tensors for 
(uncoupled) dark matter and standard model fields. The 
existence of the scalar field and its coupling change the form 
of the energy momentum tensor leading to potential changes 
in the background cosmology and structure formation. 

The coupling to a scalar field produces a direct interac- 
tion (fifth force) between dark matter particles due to the 
exchange of scalar quanta. This is best illustrated by the 
geodesic equation for dark matter particles 

where r is the position vector, t the (physical) time, $ the 
Newtonian potential and V is the spatial derivative. C v = 
dC/dip. The second term in the right hand side is the fifth 
force and only exists for coupled matter species (dark matter 
in our model). The fifth force also changes the clustering 
properties of the dark matter. 

To solve the above two equations we need to know both 
the time evolution and the spatial distribution of ip, i.e. 
we need the solutions to the scalar field equation of motion 
(EOM) 



V Vaf H j h PDM ; = 0, 



dtp 

or equivalently 
dtp 

where we have defined 
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2 THE THEORETICAL MODEL 

This section lays down the theoretical framework for investi- 
gating the effects of coupled scalar field(s) in cosmology. We 
present the relevant general field equations in § 2.1, specify 
the models analysed in this paper in § 2.2, and then briefly 
summarise the spherically symmetric solutions in § 2.3. 



2.1 Cosmology with a Coupled Scalar Field 

The equations presented in this sub-section are derived and 
discussed in (Li & Zhao 2009, 2010; Li & Barrow 2011). 
They will be used extensively in the rest of this paper and 



Veff(<p) = V(ip) + pomC(ip). 



(6) 



The background evolution of ip can be solved easily given 
the present day value of pdm since /9dm oc a . We can 
then divide ip into two parts, ip = (p + dip, where <p is the 
background value and Sip is its (not necessarily small nor 
linear) perturbation, and subtract the background part of 
the scalar field equation of motion from the full equation 
to obtain the equation of motion for Sip. In the quasi-static 
limit in which we can neglect time derivatives of Sip as com- 
pared with its spatial derivatives (which turns out to be a 
good approximation on galactic and cluster scales), we find 
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where Pdm is the background dark matter density. 

The computation of the scalar field tp from the above 
equation then completes the computation of the source term 
for the Poission equation 



V $ = - [C(tp)p OM — C(y?)pDM + 5p B 



2SV(<p)} 



(8) 



where Sps = Pb — Pb and SV(tp) = V(tp) — V(tp) are respec- 
tively the density perturbations of baryons and scalar field 
(we have neglected perturbations in the kinetic energy of the 
scalar field because it is always very small for our model). 



2.2 Specification of Model 

As mentioned above, to fully fix a model we need to specify 
the functional forms of V(tp) and C(tp). Here we will use the 
models investigated by Li & Zhao (2009, 2010); Li (2011), 
with 



C(tp) = exp^V^tp), 



and 



V{tp) = 



A 



[1 



exp l 



(9) 



(10) 



In the above A is a parameter of mass dimension four and 
is of order the present dark energy density (tp plays the role 
of dark energy in the models). 7, a are dimensionless pa- 
rameters controlling the strength of the coupling and the 
steepness of the potentials respectively. 

We shall choose a < 1 and 7 > as in Li & Zhao 
(2009, 2010), ensuring that V e // has a global minimum close 
to tp = and d 2 V e ff(tp)/dtp 2 = m 2 , at this minimum is very 
large in high density regions. There are two consequences of 
these choices of model parameters: (1) tp is trapped close to 
zero throughout cosmic history so that V(tp) ~ A behaves 
as a cosmological constant; (2) the fifth force is strongly 
suppressed in high density regions where tp acquires a large 
mass, mL 3> H 2 (H being the Hubble expansion rate), 
and thus the fifth force cannot propagate far. The sup- 
pression of the fifth force is even stonger at early times, 
thus its influence on structure formation occurs mainly at 
late times. The environment-dependent behaviour of the 
scalar field was first investigated by Khoury & Weltman 
(2004); Mota & Shaw (2007), and is often referred to as the 
'chameleon effect'. 



2.3 Solutions in Spherical Symmetric Systems 

In this subsection we summarise the solutions to the ra- 
dial profile of the scalar field tp in a spherically symmetric 
top-hat overdensity with radius Rth, and (constant) matter 
density pi n (p ou t) inside (outside) -Rth- Such a spherically 
symmetric system will be used to model dark matter halos 
later. More details concerning these solutions can be found 
in Khoury & Weltman (2004). 

If Pin = Pout, namely the mater density is the same 
everywhere, then <p will be constant across the whole space 
and its value simply minimises the effective potential V e //. 
When pi n / pout, Veff is minimised by tpi n and ip out inside 
and outside -Rth respectively, while tp will develop a non- 
trivial radial profile. 

Suppose we go towards the centre of the sphere from 



outside. If the difference between ip in and <p ou t is small, then 
tp will settle to <pi n (from tp ~ <p ou t outside) soon after we 
enter the sphere; if, on the other hand, the difference is large, 
then tp may never settle to <pi n even at the centre of the 
sphere. Khoury & Weltman (2004) give an estimate of the 
distance AR that is needed for tp to settle to <pi n from Rth: 

x/^out — s/iitp- m 



AR 

Rth 67$th 

where $th is the Newtonian potential at the surface of the 
sphere: 

k Mth k 
Z 6' 



— TTPinRTH) 



(11) 



(12) 



|7TR TH pi 



8ir Rth 

and Mth = 
sphere. Using this, Eq. 

AR lout - fin 

Rth 7\/fi;pinR, 

Khoury & Weltman (2004) present the solutions to tp 
in two regimes. In the thin-shell regime, where AR <C Rth, 
the solution is approximately 



2 • 

th 



! is the mass enclosed within the 
11) can be re-expressed as 

(13) 
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i n>2 



Vout 



AR \/^W^gTH e -m out (r-R TH ) 
-Rth t 



r e [0,R ]; 

r e [Ro,Rth].(14) 

r G [Rth, 00]; 



in which Ro G (0, Rth) and Rth — Ro <S Rth; ?riout is the 
effective mass of the scalar field outside the sphere, which is 
given by 

_ d 2 V e ff (<Pout) 



2 

m out 



dtp 2 



(15) 



In the thick-shell regime, where AR > Rth, the solution is 
approximately 




^j7$TH + |7^PinJ~ 2 , 



3r 



-m ou t(i — -Rth) 



t g [0,R TH ]; 
r G [Rth , 00 



,(16) 



Physically, if tp has developed a thin shell near the edge of 
the spherical overdensity, then from Eq. (14) we can see that 
only a fraction AR/Rth of the total mass enclosed in Rth 
contributes to the fifth force on a test particle at the edge. 
This means that the fifth force from the matter inside the 
sphere is strongly screened. In the thick-shell regime the fifth 
force is not screened. 

Note that in the thick shell regime at the edge of the 
halo we have 

|v^ 7 l[^ W ]=-2 7 ^, (17) 

which indicates that the magnitude of the fifth force is 27 2 
times that of gravity, and its effect is to rescale the Newton's 
constant by 1 + 27 2 . 



3 ANALYTICAL METHOD FOR STRUCTURE 
FORMATION 

Having reviewed the chameleon model and the solutions in 
spherically symmetric top-hat overdensities, let us now turn 
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to excursion set theory (Bond et al. 1991), which was de- 
veloped to study structure formation in cold dark matter 
scenarios. We will generalise the excursion set approach to 
the chameleon model, where the dark matter particles ex- 
perience an extra, environment-dependent, fifth force. 



3.1 Excursion Set Theory 

It is widely accepted that the large-scale structure (LSS) 
in the Universe has developed hierarchically through grav- 
itational instability. The excursion sets (regions where the 
matter density exceeds some threshold when filtered on a 
suitable scale) generally correspond to sites of formation of 
virialised structures (Schaeffer & Silk 1988; Cole & Kaiser 
1988, 1989; Efstathiou et al. 1988; Efstathiou & Rees 1988; 
Narayan & White 1987; Carlberg & Couchman 1988). 

The filtered, or smoothed, matter density perturbation 
field <5(x, R), is given by 



6(x,R) = W(\x-y\;R)5(y)d*y, 



W(k;R)5 k e Ax d% 



(18) 



where W(r;R) is a filter, or window function, with radius 
R, and W(k; R) its Fourier transform; <5(x) = p(x)/p — 1 is 
the true, unsmoothed, density perturbation field and <5k its 
Fourier transform; we will always use an overbar to denote 
background quantities. 

As usual, we assume that the initial density pertur- 
bation field <5(x) is Gaussian and specified by its power 
spectrum P{k). The root-mean-squared (rms) fluctuation of 
mass in the smoothing window is given by 



S(R) 



a 2 (R) 



(S 2 (x;R)) 



P{k)W(k-R)d 3 k. (19) 



Note that, given the power spectrum P(k), S, R and M are 
equivalent measures of the scale of a spherical perturbation 
and they will be used interchangeablly below. 

If W(k; R) is chosen to be a sharp filter in k-space, then 
the increment of 5(x; R) as R — > R — 8R or equivalently 
S —> S + 5S comes from only the extra higher-fc modes of 
the density perturbation (see Eq. (18)). The absence of cor- 
relation between these different wavenumbers means that 
the increment of <5(x; R) is independent of its previous value 
(the Markov property). It is also a Gaussian field, with zero 
mean and variance 5S. Thus, considering S as a 'time' vari- 
able, we find that <5(x; 5") can be described by a Brownian 
motion. 

The probability distribution of <5(x;_R) is a Gaussian 



P(S, S)dS ■ 



: exp 



~2S 



d8. 



(20) 



In an Einstein-de Sitter or a ACDM universe, the linear 
growth of initial density perturbations is scale-independent, 
so that 5(x) and a(R) = VS grow in the same manner, and 
as a result the density field will remain Gaussian while it is 
linear. Following the standard literature, hereafter we shall 
use <5(x; R) to denote the initial smoothed density perturba- 
tion extrapolated to the present time using linear perturba- 
tion theory, and the same for a or 5*. 

In the standard cold dark matter scenario, the initial 



smoothed densities which, extrapolated to the present time, 
equal (exceed) 5 C correspond to regions where virialised dark 
matter halos have formed today (earlier). In an Einstein-de 
Sitter universe 5 C is a constant, while in a ACDM universe it 
depends on the matter density Q m . In neither case does 5 C 
depend on the size of (or equivalently the mass enclosed in) 
the smoothed overdensity, or the environment surrounding 
the overdensity. 

As a result, to see if a spherical region with initial radius 
R has collapsed to virialised objects today or lives in some 
larger region which has collapsed earlier, we only need to see 
whether <5(x; ^ R) ^ 8 C . Put another way, the fraction of 
the total mass that is incorporated in virialised dark matter 
halos heavier than M = ^ivR^pi is just the fraction of the 
Brownian motion trajectories 5(x; S) which have crossed the 
constant barrier 5 C by the 'time' S = S(R), which is given 
by (Bond et al. 1991) 




F(M, 



where the lower limit of the integral is 



(15, 



(21) 



g+(Q) 



<5 C , because if 



d+0) 

a virialised object formed at redshift z, then its correspond- 
ing initial smoothed density linearly extrapolated to z is 5 C , 



while extrapolated to today it is 



p+(Q) 



5 C with D + (z) being 



£>+(*) 

the linear growth factor at z. In Einstein-de Sitter cosmology 
D+(z) oc (1 + z) -1 and this quantity becomes (1 + z)5 c . 

Alternatively, one can say that the fraction of the total 
mass that is incorporated in halos, the radii of which fall in 
[R,R + 5R] (or equally [S, S + 5S]) and which collapse at 
z — Zf is given by 



f(S, z f )dS = 



D+(0)5 C 



/2^sD+(z f )S 



exp 



D'i(0)5'i 



2D%{z,)S 



dS, (22) 



where f(S) the distribution of the first-crossing time of the 
Brownian motion to the barrier D+(z = 0)5 c /D+(z — Zf). 
Once this is obtained, one can compute the halo mass func- 
tion observed at Zf as 



dM 



Pm(zf) 

M 



f(S)dS. 



(23) 



dn(M) 
dM 

Other observables, such as the dark matter halo bias 
(Mo & White 1996), merger history (Lacey & Cole 1993), 
void distribution (Sheth & van de Weygaert 2004), can be 
computed with certain straightforward generalisations of the 
theory. 



3.2 Characterising the Chameleon Effect 

To incorporate the chameleon effect into the model, we need 
to have some idea about which physical quantities are most 
relevant and how they might affect the analysis. In our study 
of dark halo formation based on spherical collapse of top-hat 
overdensities, Eqs. (13, 14, 16) roughly characterise where 
the chameleon effect is strong using the following relevant 
physical quantities (in addition to the parameters a and 7 
which are fixed once a model is specified): 

(i) <fiout, the value of <p which minimises V e ff(<p) outside 
the sphere. This in turn depends on the matter density p out 
which we take approximately by smoothing the density field 
using a filter centred at our sphere with a radius £. Evidently, 
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p ou t describes the environment-dependence of the chameleon 
effect, while £ is the size of the environment, which itself is 
modelled as a spherical top-hat overdensity or underdensity. 

(ii) (fin, that minimises V e ff(ip) inside the spherical halo. 
This depends on pin, which is the density of the spherical 
halo. 

(iii) -Rth, the radius of the top- hat spherical halo. 

In summary, there are three quantities which determine the 
strength of the chameleon effect: p ou t, Pin and -Rth, of which 
the latter two characterise the spherical halo under study 
while the former represents the local environment in which 
the halo is located. 

The complexity, however, is that all of these three quan- 
tities evolve in time, and they can all be different for different 
halos. In particular, pi n , p out are the true non-linear densities 
inside and outside the halo at arbitrary redshifts z 0, while 
in the excursion set approach we are dealing with overden- 
sities which are linearly extrapolated to the present day. We 
must be able to relate the former to the latter to facilitate a 
statistical treatment based on the Gaussian distribution of 
the linear matter density perturbation field. 

3.3 Fixed-scale Environment Approximation 

In considering the linearly extrapolated matter density field, 
we must decide whether the linear evolution should be com- 
puted as in ACDM or the chameleon model? Since we as- 
sume that the chameleon model starts with the same initial 
conditions as the ACDM model, and the linear perturbation 
for the latter is much easier to compute, in what follows we 
shall always use the ACDM-linearly-extrapolated S(x;R). 

Let us consider the non-linear evolution of a smoothed 
density perturbation <5(x; R) which is surrounded by another 
top-hat sphere with ACDM-extrapolated density perturba- 
tion <5onv(x;£). It is evident that to specify the environment 
we need to know the value of f . 

There are certain guidelines in the choice of £. To rep- 
resent the local environment, £ can not be too large be- 
cause otherwise the matter density within £ would simply 
be the background value p m . £ cannot be too small either, 
because the environment should be significantly larger than 
the hosted dark matter halo to be compatible to the charac- 
teristic length scale on which the scalar field value changes 
from (pin to <£ ut . These considerations suggest that the nat- 
ural choice of £ is a few times the virial radius of the hosted 
halo. However, this means that £ is dependent on both time 
and halo size, precluding a simple analytic extension of the 
excursion set approach. 

Since we are interested in this paper in qualitative 
(rather than high precision) results, we adopt a fixed-scale 
environment approximation, in which £ is taken to be a con- 
stant. As a simple choice, we adopt £ = 8/i _1 Mpc, where 
h — Ho/ '100 km/s/Mpc and Ho is the present Hubble con- 
stant. As shown in Fig. 4 of Li & Zhao (2010), the length 
scale of the spatial variation of the scalar field value (<p ou t) is 
typically a few Mpc at late times, which is roughly the same 
as £. Such a large scale is well beyond the Compton length 
of the scalar field tp, and so the fifth force is not expected 
to play an important role. As the cosmic background expan- 
sion in the chameleon model is indistinguishable from that 
of ACDM as well, the non-linear evolution of the spherical 



overdensity enclosed by £ is well described by ACDM. This 
means that we can relate <5 cnv (x;£) to A cnv (x;£) (we shall 
use A to represent non-linear density contrasts throughout 
this paper) using the ACDM spherical collapse model, and 
then p ut = p [1 + A cnv (x;£)]. In this way, we have related 
Pout to <5 C nv(x; £). 

Assuming no shell crossing, the mass enclosed by the 
(comoving) smoothing radius R is M = ^irpR 3 . With p ou t at 
arbitrary time known, we can calculate the evolution of the 
initial density perturbation corresponding to 5(x; R) since 
(i) we know the strength of the fifth force at arbitrary time 
from Eqs. (14,16), and (ii) we can compute the collapse his- 
tory of the sphere, namely -Rth: because of mass conserva- 
tion, |7rpi n i? TH = M, giving pi n in terms of M (equivalently 
R) and -Rth , and this can be used to quantify the chameleon 
effect for the next step. 

As a result, once a top-hat overdensity S(x;R) and its 
environment 5 en v(x;£) are fixed, we can determine its col- 
lapse history. 

3.4 Spherical Collapse 

We have seen above that the spherical collapse of a 
top-hat overdensity is specified by 5(x; R) and <5 onv (x;£). 
Now we shall use these quantities to calculate the 
critical (ACDM-linearly-extrapolated) density contrast 
<5 c (x; R, Zf, <5 C nv(x; £)) that is needed for an initial overden- 
sity with radius R, residing in environment <5 env (x; £), to col- 
laspe into a virialised object at redshift zj in the chameleon 
model. In the Einstein-de Sitter and ACDM cosmologies 8 C 
does not depend on R or 5 env , but in the chameleon model 
these quantities are crucial in determining the effect of the 
fifth force. 

In the chameleon models considered here, the choice of 
parameters a and 7, as mentioned above, ensures that the 
background cosmic expansion is practically indistinguish- 
able from that of ACDM (Li & Zhao 2009). For simplicity, 
the evolution of the scale factor a(t) is specified as 

rr2 

—2 = Q m a- 3 + Q A , (24) 

H Q 

with H = a I a and the overdot denotes the (physical) time 
derivative. Throughout this paper we shall adopt Q m = 0.24, 
Qa = 0.76 and Ho = 71.9 km/s/Mpc. Also note that our 
study is limited to late times, when structure becomes non- 
linear, which is why radiation is not included in this and 
subsequent equations. 

3.4-1 Evolution of Over-densities in the ACDM Model 

Let us consider first the linear and non-linear evolution for 
an initial density perturbation in the ACDM model, which 
will be used to calculate the relation between <5 env (x; £) and 
A en v(x, a; £) (here we have written explicitly the dependence 
of A onv on time or equivalently a ox z). The convention and 
definitions here follow closely that of, e.g., Valageas (2009). 

The linear evolution of the density perturbation satis- 
fies, 

5 + 2H5 - ^Kp m S = 0. (25) 
Using equations (24) and (25), it is straightforward to show 
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that the linear growth factor D+ satisfies 



D'! + 



2 ^ £~lr, 



(N) 



D'+ - 



.(N)D+ = 0, 



(26) 



in which a prime means derivative with respect to iV = In a, 
and 



Ov(iV) 



{ 1)1 



-3N 



- 3N + Oa : 



3 -3JV 



(27) 
(28) 



are respectively the fractional densities for matter and dark 
energy at arbitrary TV. The initial conditions are given by the 
fact that, deep into the matter dominated era, D+(ai) = a^, 
and therefore D' + (at) = ai. 

To analyse non-linear spherical collapse, let us denote 
the physical radius of the considered spherical halo at time 
t by r(t), and its physical radius if it has not collapsed by 
q(t) = a(t)R (remember that R is the comoving radius of 
the filter) . Because of the spherical symmetry, it is straight- 
forward to write down the evolution equation for r(t) as 



2pa) 



(29) 



where p m = 3M/4irr s is the true matter density in the halo 
and the constant p\ is the dark energy density. Let us define 
y(t) — r(t)/q(t) and change the time variable to N. By using 
Eqs. (24, 29) and q(t) oc a(i), it can be shown that 

y' + ^#^(y- 3 -i)y = o, (30) 



y + 



2--n m (jv) 



which is clearly a non-linear equation. At very early times 
we must have y « 1 and we can write y = 1 + e with |e| <C 1. 
Substituting this into Eq. (30) to get the linearised evolution 
equation for e, and comparing with Eq. (26), we find that 
e oc D+ , in which the proportional coefficient could be found 
using mass conservation y 3 (l + Si) = 1 => e = — <5j/3 oc D+ 
(here Si is the linear density perturbation at the initial time). 
As a result, the initial conditions for y are y{a,i) = 1 — <5{/3 
and y'(a,i) = — Si/3. 

Eqs. (26, 30), associated with their corresponding initial 
conditions, completely determine the necessary dynamics in 
the ACDM model. 



3.4- 2 Evolution of Overdensities in the Chameleon Model 

With the preliminaries given above, we can now consider 
spherical collapse in the chameleon model. 

From the discussions in § 2.3 and results of Li & Zhao 
(2009), we know that the fifth force acts as if it renormalises 
the Newton's constant by (I + 27 2 ) if the chameleon effect is 
weak (i.e., in the thick-shell regime); on the other hand, it is 
strongly suppressed in the thin-shell regime. In particular, 
comparison of Eqs. (14, 16) shows that the two regimes give 
the same exterior solution when AR/Rth = 1/3. Therefore, 
we propose to approximately take account of the effect of the 
fifth force as if it effectively rescales the Newton's constant 
by I + 27 2 min {3AJ?/_Rth, 1}- This is certainly not expected 
to be very accurate, but our aim here is to present a method 
which captures the essential features of the environment de- 
pendence. 

Because we do not need the the linear perturbation evo- 
lution in the chameleon model, we shall go to the spherical 



collapse directly. According to the above approximation, the 
equation of motion of a spherical shell at the edge of the top- 
hat overdensity is 



1 



- = g K PA - "Kpm 1 + 27 min {"^^"' 1 } ' ( 31 ) 

where we have neglected the perturbation in the energy den- 
sity of the scalar field and its kinetic energy, which are neg- 
ligible (Li & Zhao 2009). Note that this means that the en- 
ergy density of the scalar field is the same as that of the 
vacuum energy in ACDM model. 

The scalar field value which minimises the effective po- 
tential V e ff(<p) is given by (Li & Zhao 2009) 

a Vo 

7 Pr, 

where p m is the local matter density. Substituting this into 
Eq. (11), we find that, at time a, 

3AR 1 a Q A 



r- a Vq 

VKip W , 



(32) 



Rth 



(H0R) 2 7 2 



9 2 

- - 1. 



Vcnv 



1 + S C 



1 + 5, 



y h a , (33) 



in which j/h is the y for the considered halo and y cn v that for 
the environmental spherical overdensity smoothed at radius 
£• 5 e nv,i and Si are respectively the initial values for 5 en v(x; £) 
and S(x; R) and 

D + (z = Zi 



D+{z = 0) 
D + (z = Zi) 



6(x;R). 



(34) 



D+{z = Q) 

In the derivation of Eq. (33) we have used the approximation 
that masses are conserved within the top-hat overdensities 
with radii R and £. Note that because of the unit convention 
c — 1 the quantity HqR is dimensionless. Eq. (33) shows that 
the effects of the fifth force will be more suppressed by: 

(i) increasing 7 and decreasing a, both making the scalar 
field heavier and unable to propagate far; 

(ii) increasing Q m , meaning that the matter density is 
higher in the Universe, again making the scalar field heavier; 

(iii) increasing environmental density 5 en v(x; £), therefore 
making the term in the brackets smaller; 

(iv) considering earlier times, where a is smaller, because 
the overall matter density is higher then, and 

(v) considering bigger halos (larger R), which are more 
efficient in screening the fifth force. 

From the earlier discussion, y env is governed by Eq. (30), 
and now we need to find an evolution equation for yh as well. 
This can be obtained from Eq. (31) following the derivation 
of Eq. (30). The result is 

3, 



" 1 
Vh + 


2 - 




[N) 



rOn(iV) 



Vh 



1+27 Hiw 1 ] 



(35) 



where 3A_R/_Rth is given by Eq. (33). Because at very early 
times the chameleon effect is very strong, the initial condi- 
tions of this equation can be chosen exactly as in the ACDM 
model. Eqs. (26, 30, 35), together with Eqs. (33, 34) form 
a closed system for our chameleon model. These completely 
fix the evolution of a spherical overdensity <5(x; R) residing 
in the environment <5 cnv (x;£). Note that Eq. (26) only needs 
to be solved once. 
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Figure 1. The critical (ACDM-lincarly-cxtrapolated) density 
perturbation S c for the given spherical overdensity with mass M 
to collapse at Zf = 0. Shown are <5 C as functions of M for halos 
residing in different environments, with (solid curves from top to 
bottom) S env = 1.6, 1.2, 0.8, 0.4, 0, -0.4, -0.8, -1.2. For compari- 
son, the constant 6 C ~ 1.676 for the ACDM model is overplottcd 
as the dashed horizontal line. 
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Figure 2. The critical (ACDM-linearly-extrapolated) density 
perturbation S c for the given spherical overdensity residing in en- 
vironment <5 e nv to collapse at Zf = 0. Shown are 5 C as functions 
of <5env for halos with different masses (solid curves from top to 
bottom) M = 10 15 , 10 14 , 10 13 , 10 12 , 10 11 h- 1 M SUB . For compari- 
son, the constant 8 C ~ 1.676 for the ACDM model is overplottcd 
as the dashed horizontal line. 



3.4-3 Numerical Examples 

To get an idea about how the environment-dependent fifth 
force changes spherical collapse in the chameleon model, we 
present some numerical examples in this section. 

As we have discussed above, the critical density (linearly 
extrapolated to today using ACDM model) which is needed 
for a spherical overdensity to collapse at redshift Zf depends 
on R (the spherical overdensity's own property) and 5 env (its 
environment): <5 c (x) = 5 c (x; R, 6 eav , Zf) = <5 c (x; M, 5 env , Zf) 
where we have used M « |7rp m o_R 3 . 

Fig. 1 shows S c as a function of halo mass M for different 
values of <5 cnv and Zf = 0. We have considered halos in eight 



Figure 3. Same as Fig. 1 but for spherical overdensities which 
collapse at zj = 1. In this case <5 C fs 2.657 for the ACDM model 
(the horizontal dashed line). 



different environments with 8 cnv ranging between 1.6 (very 
dense environment) and —1.2 (very empty envrionment). As 
can be seen there, the fifth force lowers S c compared with the 
ACDM result (dashed line) , which is as expected because it 
makes collapse easier. Note that 

(i) Unlike in ACDM, in the chameleon model S c is mass 
and therefore scale dependent, a point which we will return 
to later. 

(ii) For a given <5 en v, S c is closer to the ACDM result for 
more massive halos because these halos are more efficient 
in screening the fifth force [see also Eq. (13)]. Note however 
that S c will never exceed the corresponding value in ACDM 
model because the fifth force always helps rather than pre- 
vents the collapse. 

(iii) For a given halo mass M , 8 C is closer to the ACDM 
prediction in denser environments, where the chameleon ef- 
fect is stronger. 

These can also be seen in Fig. 2, which shows 8 C as a function 
of S cnv for different halo masses. 

Fig. 3 shows the same results as Fig. 1, but for the ha- 
los which collapse at Zf = 1. This shows similar qualitative 
behaviour as does the Zf = case, but the relative dif- 
ference between the collapsing threshold <5 C and its ACDM 
result is generally smaller because by z = 1 the fifth force 
is strongly suppressed by the chameleon mechanism in most 
environments and because the halos which form at z = 
experience the fifth force for longer. 

The simplified computation described in this section can 
capture the essential effects of the chameleon fifth force. We 
will use it as an ingredient of the extended excursion set 
model to be introduced below. 



3-4-4 Notes on the Approximations Used 

As mentioned earlier, the purpose of this work is to intro- 
duce a conceptually simple, largely analytic, method of in- 
corporating environment dependence in the study of struc- 
ture formation that is adequate for parameter exploration. 
Consequently we have used a number of approximations to 
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simplify the calculation. Here we briefly summarise these ap- 
proximations and discuss how they can be improved using 
numerical methods: 

(i) The computation of the scalar field profile <p(r) in the 
spherical halo: in this work we have adopted the analytical 
approximations given in Khoury & Weltman (2004), which 
could be improved by solving the scalar field EOM explicitly 
using numerical methods. 

(ii) The detailed shape of the spherical halo: because of 
the environment dependence of the fifth force, shells at dif- 
ferent radii of the halo will travel at different speeds, result- 
ing in a modification to the top-hat shape of the halo. In this 
work we have assumed a constant overdensity for the halo, 
which is only an approximation. In general we expect that 
matter will accumulate (slightly) towards the edge of the 
halo. This effect can be computed accurately once <p(r), or 
equivalently the profile of the fifth force is known precisely 
(see Martino et al. (2009) for an example). 

We will leave these improvements to future work. 



3.5 Generalised Excursion Set Method for the 
Chameleon Model 

We have seen above that the excursion set prediction of the 
halo mass function (based on the spherical collapse model in 
the ACDM cosmology) is closely related to the first crossing 
distribution of a flat barrier by a Brownian random walk 
that starts from zero. In the chameleon model two factors 
lead to a more complicated problem. 

(i) The barrier that is to be crossed by the Brownian mo- 
tion is no longer flat, but rather depends on the mass scale 
M (c.f. Figs. 1 and 3), or equivalently R or S(R). 

(ii) The barrier is also affected by the environment sur- 
rounding the collapsing halo (c.f. Fig. 2), and so we need to 
know the probability distribution of its environment (S env ) 
as well. 

These complications are the subject of this section. 



3.5.1 Unconditional First Crossing of a Moving Barrier 

The distribution of the first crossing of a general barrier 
by a Brownian motion has no closed-form analytical solu- 
tions except for some simple barriers, e.g., flat (Bond et al. 
1991) and linear (Sheth 1998; Sheth & Tormen 2002). Un- 
fortunately neither of these is a good approximation to our 
general barrier (cf. Fig. 1). As a result, we shall follow 
Zhang & Hui (2006) and numerically compute this distribu- 
tion. We shall briefly review their method for completeness. 

Denote the unconditional probability that a Brownian 
motion starting off at zero hits the barrier 5 C (S) for the first 
time in [S, S + dS] by f(S)dS. Then, f(S), the probability 
density, satisfies the following integral equation 



f(S) = g(S) + / dS'f(S')h(S,S'), 
Jo 

in which 



(36) 



g(S) 



8 C _ dS c 
~S~~dS 



P(Sc,S).. 



= 1CT', 7 = 0.5, z, = 0.0 




Figure 4. (Colour online) The moving barriers 8 C (S) for different 
values of 5 env as indicated beside the solid curves. The dashed line 
is the constant 8 C for spherical collapse in the ACDM model. The 
vertical dotted line represents S = = cr| = 0.64, which gives 
the length scale used to define environment. Also plotted is the 
trajectory of a Brownian random walk which starts at (Sc, S cnv = 
1.0) (the triangle). Note that the first crossing happens earlier in 
the chameleon models because the barrier is lower. 



h(S,S') = 



,dS c 



dS S- 



5' 



P(S C -S' C ,3-S'), (37) 



where for brevity we have suppressed the S'-dependence of 
S C (S) and used S' c = 5 C (S'); P(S, S) is given in Eq. (20). This 
equation could be solved numerically on an equally-spaced 
mesh on S: Si = iAS with i = 0, 1, • • • , N and AS = S/N. 
The solution is (Zhang & Hui 2006) 



fa = 9o = 0, 

/i = (l-A 1A )- 1 g 1 , 



(38) 



= (1-Ai, 



3=1 



where we have used ft = f(Si) and similarly for gi to lighten 
the notation, and defined 



AS 



(Si, 



AS \ 
2 ) 



(39) 



We have checked that this method agrees accurately with 
the analytic solution for the flat-barrier crossing problem. 



3.5.2 Conditional First Crossing of a Moving Barrier 

The unconditional first crossing distribution, which relates 
directly to the halo mass function in the ACDM model, is 
not particularly useful in the chameleon model. This is be- 
cause spherical overdensities in different environments will 
follow different evolution paths. If it is in the environment 
specified by (<5 onv , St), then (<5 env , S%) should be the starting 
point of the Brownian motion trajectory. In other words, we 
actually require the first crossing distribution conditional on 
the trajectory passing <5 onv at S = S( . Note that in a broader 
sense the unconditional distribution is a conditional one with 

(5env,S ( ) = (0,0). 

Evidently, <5 cnv has its own distribution: very dense and 
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Figure 5. (Colour online) The condition first-crossing distribu- 
tion for Brownian random walks starting off at (S^,8 en \r) in the 
chameleon (the solid curves) and ACDM models (dashed curves). 
This depends sensitively on the values of <5 e nv (indicated above 
the curves) as explained in the text. Some physical parameters 
arc also shown. 
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Figure 6. Upper panel: the mass functions for the chameleon 
(solid curve) and ACDM (dashed curve) models. Some physical 
parameters are shown in the figure, and others include Q m = 
0.24. Lower panel: the fractional differences between the two mass 
functions (solid curve); the dashed line is identically zero and is 
shown as a reference. 



very empty environments are both quite rare. To quantify 
this distribution, we need to first define the environment, or 
equally its smoothing scale £, which has been chosen to be 
8/i _1 Mpc above. 

The problem then reduces to the calculation of the first 
crossing probability conditional on the Brownian motion tra- 
jectory passing S cnv at S$ = erf: f(S c (S, S cnv ), S \ S crLV ,S ( ), 
where we have written explicitly the 5 en v-dependence of 8 C . 
The numerical algorithm to calculate the conditional first 
crossing probability is a simple generalisation of the one used 
above to compute the unconditional first crossing probabil- 
ity (Parfrey et al. 2011) and is not presented in detail here. 

Fig. 4 shows the moving barrier 5 C (S) as a function of 
S for different values of <5 e nv As an illustration, we have also 
shown a Brownian motion trajectory which passes <5 cnv = 1.0 
at = erf = 0.64 (the triangle). Clearly, the larger the 
value of S cnv , the more likely the Brownian motion will hit 
the barrier at smaller S. This is what we see in Fig. 5, which 
shows the conditional distribution f(5 c (S, <5 cnv ), S | S cnv , £{) 
for different values of <5 cnv . 

For comparison we also show the corresponding results 
for the ACDM model using the dashed curves in Fig. 5. 
Note that the solid curves are always higher than the dashed 
ones for smaller S and lower for bigger 5*. This is because 
in the chameleon model the barrier is generally lower and 
the Brownian motion is likely to cross it for the first time at 
smaller S. 



3.5.3 Integrating over the Environment Distribution 

To get the final first crossing distribution of the moving bar- 
rier, we must integrate over all environments. The distribu- 
tion of <5 e nv, denoted as q(8 cnv , 8 BC , £{), in which 5 SC is the 
critical overdensity for the spherical collapse in the ACDM 



model 1 , is simply the probability that the Brownian motion 
passes S cnv at S% and never exceeds S sc for S < S% (because 
otherwise the environment itself has collapsed already) . This 
has been derived by Bond et al. (1991): 



i 



exp 



2 

env 



exp 



(8 C 



2<5 SC ) 



2^ 



(40) 



for 8 env ^ S BC and otherwise. 

Then the environment-averaged first crossing distribu- 
tion will be 



/ave {S) 



q x f(8 c (S, 8 

env ) i $ | ^env j S f )d(W (41) 



In the special case where the barrier is flat, 8 C (S, S en v) = S sc , 
f(8 c (S,S anv ),S | S e nv,S^) is known analytically as 



/ = 



2tt (5 - S e ) 8/2 



exp 



(^sc ^env ) 

' 2(5*-5 e ) 



(42) 



and the integration in Eq. (41) can be performed exactly to 
obtain 



fave{S) 



■ exp 



Si 
2S 



(43) 



which is just the unconditional first crossing distribution for 
a constant barrier 8 SC at S. This is as expected, because the 
collapse does not depend on the environment. 

In general cases with environment-dependent collapse, 
/ave(S) must be computed numerically. Indeed, in Eq. (41) 
both <?(<5env, Ssc, S^) and f(S c (S, cWv), S | e> cn v, S() differ from 
the flat-barrier case. The distribution / has been discussed 
above (cf. Fig. 5). The distribution q should, in principle, 



1 Remember again that the evolution of the environment is as- 
sumed to be governed by the ACDM model. 
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be calculated for the chameleon model numerically, but we 
choose to use the ACDM result Eq. (40) for the following 
reasons: recall that g(<5 cnv , 8 SC , S$) is the probability that the 
Brownian motion starts off at the origin, never hits the con- 
stant barrier S sc before and goes through 8 en v at S%. To 
estimate its difference from the true value in the chameleon 
model, we replace the S sc in Eq. (40) with 5 C (S^) for the 
<5 C nv values in the figures, and find that the change of q is at 
percent and subpercent level 2 , which is not surprising given 
that S C (S ^ S$) is very close to S sc (cf. Fig. 4) (the approx- 
imation will be even better for higher redshift cf. Fig. 3). A 
better approximation would be to assume S C (S) m 8 — j3S 
where /3 is some constant, but here we do not see the neces- 
sity for doing this. 

Using Eq. (40), we perform the integral in Eq. (41) us- 
ing Gaussian quadrature. We checked the accuracy of this 
method by applying it to the flat-barrier case and find that 
the agreement with the exact solution is excellent. The halo 
mass function is related to the averaged first-crossing distri- 
bution /avc(S') by 



dS 
dM 



dM, 



and we have plotted in Fig. 6 the function dn(M, z = 
0) /dM for both the chameleon (solid curves) and the ACDM 
(dashed curve) models. As can be seen clearly, the fifth 
force results in more massive halos than in the ACDM 
model, but in compensation there are fewer low mass ha- 
los (M < 10 12 1i~ 1 A1q) in the chameleon model. 

To see the difference more clearly, we have also plotted 
the fractional difference between the chameleon and ACDM 
mass functions in the lower panel of Fig. 6. This shows that 
the increase of n(M) is largest for halos in the mass range 

10 13 < M I ^MQj < 10 14 . For high mass halos the 

fifth force is strongly suppressed and its effect on the mass 
function is smaller, as expected. Because a larger fraction 
of the total mass has been assembled in high mass halos, 
fewer small isolated halos survive the merger and accretion 
process 3 . 

Note that the effects of the fifth force are suppressed for 
high mass halos, not only because the halos are efficient at 
screening that force themselves, but also because they are 
more likely to reside in dense environments. More explicitly, 
the probability distribution of S anv at S — S$, given that the 
Brownian motion goes through S ~ S C (S) at S (where it is 
about to cross the barrier), is 



p(Scnv | S,S) 



: exp 



s e 



2^(5- 



Se 



For high mass halos, S is close to S$ and this distribution 
strongly peaks at <5 cnv ~ S^S/S ~ S. The results, of course, 



2 Noe that this is an upper bound of the error of using Eq. (40), 
because the barrier does not stay at 8 C (S^) for all S £ [0, St] but 
rather decreases from <5 SC at S = to it at S = , 

3 We want to emphasise that the result here is not directly com- 
parable with that obtained in Li & Zhao (2010) using full TV-body 
simulations, because there subhalos are also counted. 



Table 1. The parameters a and fi for the 9 chameleon models 
studied in this section. The ACDM paradigm corresponds to a = 
= 0. We assume O m = 0.24 and U A = 0.76. 

lg(ct) -7 -7 -7 -6 -6 -6 -5 -5 -5 
P 1/3 1/2 1 1/3 1/2 1 1/3 1/2 1/1 



are consistent with our intuitive understanding about the 
chameleon effect. 



4 APPLICATIONS 

The key new concept in our extended excursion set model is 
the specification of the environment in terms of two param- 
eters (S^,5 env ): the environment determines how spherical 
collapse in the chameleon model is modified compared to 
ACDM, and also means that we have to use conditional dis- 
tribution of the first crossing rather than the unconditional 
distriubtion, as in the conventional excursion set approach, 
to compute observables such as the mass function of non- 
linear structures. 

The use of the conditional first-crossing distribution is 
not new. Mo & White (1996), for example, used it to study 
the bias between the halo number density and the dark mat- 
ter density fields in the ACDM model. In the chameleon 
model, both bias and the mass function must be computed 
using the conditional first-crossing distribution. In fact, com- 
putation of the mass function is more complicated since we 
need to average over the probability distribution of environ- 
ments. 

The methods introduced here can be used to study, for 
example, the formation redshift of halos zj and their depen- 
dence on the parameters a and 7 describing the chameleon 
mechanism, Zf. The most difficult step in such a calculation 
is the computation of the moving and environment depen- 
dent barrier S C (S, 8 cnv ). Nevertheless, the computations are 
very much faster that A-body simulations and so we can 
explore large regions of parameter space rapidly. 

Let us consider the mass functions of the chameleon 
models with model parameters as specified in Table 4. Fig. 7 
shows the critical density for a spherical overdensity to col- 
lapse at Zf — as a function of the enclosed mass M and 
environment 5 env - As expected, the collapse threshold is 
lower in all chameleon models because the fifth force, how- 
ever weak, is always attractive and boosts the collapse. For 
smaller a (left column, a = 10~ 7 ), the difference from the 
ACDM prediction is small, especially for the largest over- 
densities because the fifth force is more strongly suppressed 
in these systems, as discussed in the previous section. On the 
other hand, for large a (middle and right columns, a = 10 -6 
and 10 ), the deviation from ACDM is much larger, even 
for the largest overdensities. Increasing 7 will strengthen the 
fifth force and therefore also lower the collapse threshold. 

Fig. 8 is equivalent to Fig. 7, but for spherical overden- 
sities which collapse at Zf = 1. Because the matter density 
is higher at higher redshift, the fifth force is more strongly 
suppressed and hence the deviation from ACDM is smaller. 

Finally, we have plotted the effect of a chameleon-type 
fifth force on the dark matter halo mass functions in Fig. 9. 
For clarity Fig. 9 shows the fractional change of the quantity 
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Figure 7. The critical ACDM-linearly-evolved overdcnsity for spherical collapse at zt = as a function of the mass enclosed and the 
environment <5 en v The physical parameters 0,7 are indicated beside each panel. In each panel the solid curves from top to bottom are 
respectively for 5 cnv = 1.6, 1.2, 0.8, 0.4, 0.0, -0.4, -0.8 and -1.2. The dashed line is the result for ACDM model. 




Figure 8. The same as Fig. 8, but for collapsing redshift Zj = 1. In each panel the solid curves from top to bottom are respectively for 
<5 cnv = 2.4, 1.8, 1.2, 0.6, 0.0, -0.6, -1.2 and -1.8. 
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Figure 9. (Colour online) The fractional difference of dn/dM between the chameleon (the scalar field parameters a and 7 are indicated 
beside each panel) and ACDM models, at three redshifts (black curves), 1 (green curves) and 2 (red curves). The result for ACDM is 
plotted as the dashed line for reference. 



dn/dM, where n(M) is the halo mass function, with respect 
to the ACDM prediction, at three redshifts z = 0, 1, 2 respec- 
tively. For a = 10~ 7 , the fifth force is strongly suppressed 
and the fractional change of dn/dM is less than 10%, even 
for 7 = 1 and z — 0. 

For a = 1CP 6 , the fifth force is less suppressed and 
the fractional change of dn/dM at 2 = could be up to 
~ 15% (for 7 = 1/3) or even ~ 40% (for 7=1), showing 
interesting and potentially observable effects. The deviation 
at early times is mainly restricted to lower mass halos. At 
later times, massive halos also start to feel the fifth force, and 
a deviation is seen at higher halo masses. With a — 10~ 5 , 
the qualitative features mentioned above all remain, but the 
deviation from ACDM is much stronger, up to ~ 60% for 
7 = 1/3 and more than 100% for 7 = 1 at z — 0. 

Figs. 7-9 show that some choices of parameters can 
lead to large deviations in the abundances of non-linear 
objects compared to the ACDM model. The tightest con- 
straints on the parameters a and /3 would probably come 
from number counts and number densities of well char- 
acterised galaxy cluster samples intermediate redshifts, 
z < 1. Such samples are becoming available from com- 
bined Sunayev-Zeldovich/X-ray measurements from Planck 
(Planck collaboration 2011) and SPT (Carlstrom et al. 
2011). A detailed comparison of observations with our model 
is beyond the scope of this paper. 

One question one might ask is if the chameleon model 



could be used to produce more halos at very early times, 
say z > 6, which might ease problems in reionizing the 
intergalactic medium at early times. It has been argued 
(Hellwing et al. (2010)) that a fifth force might boost hier- 
archical structure formation leading to enhanced production 
of UV photons at early times. Unfortunately, the chameleon- 
type fifth force is strongly suppressed at earlier times. Fig. 9 
shows this up to z = 2, and for z > 6 the deviation from 
ACDM is even smaller. 



5 SUMMARY AND CONCLUSIONS 

To summarise, in this paper we have presented an extension 
of the standard excursion set theory so that it can be used to 
study structure-formation scenarios which are environmen- 
tally dependent. Our method separates the calculation into 
two steps: 

(i) compute the collapse of the spherical overdensity in a 
given environment; 

(ii) compute the probability that the spherical overden- 
sity is located in the specified environment, and average over 
the distribution of environments. 

For (i) we have proposed a simplified model, which is a gen- 
eralisation of the usual spherical collapse model to the case 
in which the overdensity evolves inside an evolving envi- 
ronment. For (ii) we have derived an approximation to the 
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environment distribution, and shown how to compute the 
averaged first-crossing distribution, which is closely related 
to the halo mass functions. 

As a working example, we have applied the method to 
the chameleon model. Our numerical results agree with how 
we expect the chameleon effect to behave as a function of 
the model parameters. We have concentrated here on the col- 
lapse redshift and mass functions of virialized objects. These 
predictions could be used in conjunction with forthcoming 
data to set constraints on the model parameters. 

As in the standard excursion set theory for the cold 
dark matter model, it is straightforward to generalise the 
excursion set approach to compute other observables, for 
example the formation of voids (Sheth & van de Weygaert 
2004) and the merger history of halos (Lacey & Cole 1993). 
In conjunction with the halo model, our method can be used 
to predict the non-linear matter power spectrum as well. 

Structure formation scenarios with strong environment 
dependence have become more and more popular recently. 
Newtonian gravity has been tested to high precision in our 
local environment, but deviations may be significant on cos- 
mological scales, perhaps as a result of a chameleon-like 
mechanism. The methods presented here offer a faster alter- 
native to iV-body simulations, enabling a much wider range 
of models to be confronted with observations. Meanwhile, 
the analytic formulae presented here enable a clear track of 
the underlying physics, such as in which ways the chameleon 
effect modifies the structure formation. We hope that these 
methods will contribute towards a better understanding of 
the nature of dark energy. 
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The MNRAS KTjhX 2e class file is implemented by placing the command 
\documentclass [options] {mn2e} 
at the start of the document. 

The various option commands listed in the Monthly Notices DTJrX style guide for authors can still be used. The following 
additional options exist. 

(i) useAMS - this enables the production of upright Greek characters n, \i and 9 ($\upi$, $\umu$ and $\upartial$) and 
slanted ^ and ^ ($\leq$ and $\geq$). Characters 7t, [i and 3 appear upright only on systems that have the Euler roman 
fonts (eurmxx); characters ^ and ^ appear slanted only on systems that have the AMS series A fonts (msamxx). On systems 
that do not have these fonts, the standard forms of the characters appear in the printout; however, they should be correct in 
the final typeset paper if the correct DTpX commands have been used. 

(ii) usedcolumn - this uses the package file dcolumn.sty to define two new types of column alignment for use in tables. If 
the usedcolumn option has been specified then, within a table, d{_x} can be used to produce a 'flush right' decimal-aligned 
column with x decimal places and . can be used to provide a decimal- aligned column centred on the decimal point (the 
number of decimal places does not need to be specified for this column type). Note that the standard DTpX 'tools' packages 
dcolumn.sty and array. sty are required in order to use the usedcolumn option. 

(iii) usenatbib - this uses Patrick Daly's natbib.sty package for cross-referencing. If the usenatbib option is specified, 
citations in the text should be in one of the following forms (or one of the additional forms documented within natbib . sty 
itself). 

• \citet{fce?/} produces text citations, e.g. Jones et al. (1990), 

• \citep{fce?/} produces citations in parentheses, e.g. (Jones et al. 1990), 

• \citealt{fce?/} produces citations with no parentheses, e.g. Jones et al. 1990. 

For three-author papers, a full author list can be forced by putting a * just before the {. To add notes within the citation, use 
the form \citep \_pre_reference_texf] lpost_reference_text] {key} (note that either of pre_reference_text and postjreJerence_text 
can be blank). 

Items in the reference list must be of the form 
\bibitem [\protect\citeauthoryear{a«t/ior_names}-{year}-] {.key} Text of reference ... 
for one-, two- and multi-author papers, or 

\bibitem [\pTotect\cite&uthorye&r{three_author_namesy{.first_author_etaiy{.yeary~\ {key} Text of reference ... 
for three-author papers. 

Note that Patrick Daly's package natbib.sty is required in order to use the usenatbib option. 

We recommend that authors use natbib . sty as their standard cross-referencing package, because of the flexibility in citation 
style that it provides. 

(iv) usegraphicx - this loads the graphics package graphicx . sty, which authors can use to include figures in their papers. 
Note that the standard I^TpX graphics package graphicx. sty is required in order to use the usegraphicx option. 

Other style files, or packages providing features such as graphics support, can be used in conjunction with mn.cls. To do 
this, the command \usepackage-Cpacfcage_name} is used. 

An additional feature of the class file is that footnote symbols are now in the correct journal style (symbols for title page 
footnotes, superscript arabic numbers for text footnotes). 

A BibTeX style file, mn2e.bst, is now available for authors who wish to use BiBTeX. It is recommended that this should 
be used in conjunction with natbib.sty. 

For general instructions on preparing Monthly Notices papers in I^TffX, please refer to the MNRAS DTpXstyle guide for 
authors, available on the CTAN sites in subdirectory /tex-archive/macros/latex209/mnras. 



